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Abstract — In the design of quantum computing devices of 
the future the basic element is the qubit. It is a two-level 
quantum system which may describe population transfer from 
one steady-state to another controlled by a coherent laser field. 
A four-dimensional real-variable differential equation model is 
constructed from the complex-valued two-level model describ- 
ing the wave function of the system. The state transition matrix 
of the model is constructed via the Wei-Norman technique and 
Lie algebraic methodology. The idea of parametrization using 
flatness-based control, is applied to construct feasible input- 
output pairs of the model. This input drives the state of the 
system from the given initial state to the given final state in 
a finite time producing the corresponding output of the pair. 
The population transfer is obtained by nullifying part of the 
state vector via careful selection of the parameter functions. A 
preliminary simulation study completes the paper. 

I. INTRODUCTION 

In quantum mechanical framework deterministic bits "1" 
and "0" are substituted by the qubit [22]. The qubit is a 
composition of the pure states "1" and "0". This composition 
means that the actual state of the qubit is not exactly "1" or 
"0" but a combination of these. In the measurement, however, 
the outcome is always one of the two possibilities "1" or "0". 
The qubit can be represented as a point on the surface of a 
sphere, so-called Bloch sphere, see e.g. [22]. If one wants 
to save information into a qubit, then the key peoblem is to 
drive the qubit from one state to another. Then one arrives 
at the description of the qubit as a dynamic differential 
equation system, the controls of which are the parameters 
of the driving laser field. 

In quantum computation the qubit forms a basic element 
for building up multi-qubit computing elements of future 
quantum computers, see [16]. Then a key problem is to drive 
the qubit from one stable level to another. 

Molecular excitation, i.e. driving of an ensemble of 
molecules from one locally stable steady state to another is 
one alternative for a qubit structure. This type of systems are 
controlled by using coherent light. Based on laser technology 
shorter and shorter coherent pulses can be generated for 
controlling molecular excitation, see [16]— [19]. The goal 
is to direct molecular reactions towards unprobable but 
desirable direction [3]-[5]. Then nonlinear and more and 
more sophisticated control methods are needed for properly 
designing durations and forms of the control pulses. In 
classical N-level problems the system to be controlled can 
be modelled by using ordinary 27V-dimensional differential 
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equation systems. Due to femto- and picosecond scale pulses 
feedback is not in general applicable in the control design 
for these systems. Flatness-based control, see [6]-[9] & [11], 
is then an ideal methodology for open-loop design relevant 
in quantum control problems. 

Due to the fact that the bilinear quantum control systems 
are not controllable in the whole Euclidian space E 2Af the 
methodology applied here is called parametrization. 

This two-level quantum control problem and some related 
studies have been carried out by several authors earlier, too, 
see [1], [10], & [20]. Especially, in [20] a very similar 
approach as ours is used. 

However, we start from the basic definition of differential 
flatness. The system 

fir 

— = f(x,u); x(t)eR n , u(t)€R m (1) 

is called differentially flat if there exists algebraic functions 
([6]) A, B, C, and finite integers a, (3, and 7 such that for any 
pair (x, u) of inputs and controls, satisfying the dynamics Q, 
there exists a function z, called a flat (or linearizing) output, 
such that the following equations are satisfied 

x(t) = A(z,z, ...,^ (a) ) 

u(t)=B(z,z, ...,z^) (2) 

z(t) = C(x, it, ii, . . . , u^). 

The actual output y, which is not present in the definition of 
flatness, may have the dependence 

y(t) = h(x(t),u{t)) 

for some given output function h. In parametrization pro- 
cedure, due to uncontrollability, the last equation in ([2]) for 
z(t) is neither constructed nor applied. 

From the standard finite-state Schrodinger equation of two 
energy levels a four-dimensional real-variable differential 
equation model is obtained. The Wei-Norman technique is 
used in the construction according to [21], The exponential 
representation of the transition matrix of the system includes 
three base functions, two of which serve as the parameter 
functions. In this framework the initial and final states can 
be defined corresponding to the two levels of the original 
system model. Then parametrization design is applied for 
explicitly calculating the parameter functions, which in turn 
give the desired input-output pairs. 

II. SYSTEM MODELS 

Population transfer in a two-level quantum system, see 
[4], can be described by the time-dependent Schrodinger 



i d *t = H{t)$, H(t) 



(3) 



equation, i.e. by the dynamics 

E x O(i) 

n*(t) e 2 

where the modified Planck's constant H = ^ has been scaled 

^ C 2 has 



to h — 1, and i = v — 1- The wavefunction ip : R — 
the probabilistic interpretation, in the sense that 

||#)|| 2 H^(*)l 2 + hMt)| 2 = i,VieK, (4) 

where -0 = (ipi,ip2)- The control is given by ft : R — > C, 
and f2* is the complex conjugate of Q. and i? 2 are the 
energy levels. The unitary transformation ip i-> V an d i-> u 
by 



C/(i) 
«(*) 

transforms Q to 

.# 





i£ 2 * 
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- l (B»-BO*n( t ) 
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u*(t) 

The componentwise representation 

1 



o 
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converts <|8j to the dynamics 

ipi = -info, 

tp 2 = -iu* Ipi. 
By using the real-valued decompositions 

4>i = Xi + ix 2 
i^2 = x 3 + ix 4 
u = U\ + iu 2 
one obtains a state-variable representation 
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or in another form 

dir 
dt 
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(16) 
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The constraint Q is converted into the form 
E*2 



1. 



(18) 



fc=i 



Remark 1: The matrices F\ and F 2 together with their 
Lie product 2F 3 = [Fi,F 2 ] = FiF 2 - F 2 Fi form a Lie 
algebra. This can be used as a basis for differential geometric 



considerations of the control system ( 13 1. However, the 
elementary approach applied in this paper is sufficient for 
our parametrization purposes. 

III. WEI-NORMAN REPRESENTATION 

The Lie algebra of the matrices F\, F 2 , and F 3 is three- 
dimensional with the relations 



[F U F 2 ]=2F 3 , 
[F 2 ,F 3 ]=2F 1 , 
[F 3 ,F 1 ]=2F 2 . 



F 3 
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(19) 
(20) 
(21) 

(22) 



Due to the linear structure of the system model ( 14 1 with 
respect to the state x, the state transition matrix of the 
system, denoted by $, and which relates the values of the 
state according to 

x{t) =$(t,0)x(0) (23) 



can be written as a product of exponentials 

*(i,0) 



q9iF\ gttefb e 93^3 



(24) 



where the exponentials are defined by the absolutely con- 
verging infinite series 



k=0 



o k F k 



% = 1,2,3. 



(25) 



The state transition matrix satisafies the following initial- 
value problem (IVP1) 

9 -$(t,0) = F(t)*(t,0); $(0,0)=/, (26) 



dt 



F(t) = ui(t)F x + u 2 {t)F 2 + • F 3 



(27) 



The technique we are using is nowadays called Wei-Norman 
technique according to the paper of Wei and Norman [21]. 



Substitution of the (|24|i to the IVP1 gives 
d 



dt 



$ = <h Fi $ + g 2 e 9lFl F 2 e~ 9lFl $ 



(28) 



By using (several times) the Campbell-Baker-Hausdorff 
formula for square matrices A and B of the same dimension 

= B + [A, B] + [A, [A, B}]/2\ (29) 
+ [A,[A,[A,B]}}/3\ + --- 



e A Bc- A 



in the equation ( 28 1 it can be represented in the form 



d 

— $ 
dt 


= [hi^F^ f 2 it)F 2 + h(t)Fz] $ 


(30) 


hit) 


= 9i + 93 sin(2g 2 ) 


(31) 


hit) 


= g 2 cos(25i) - g 3 cos(25 2 ) sin(2fln) 


(32) 


hit) 


= g 2 sin(2fln) + g 3 cos(2gi 2 ) cos(2gi) 


(33) 



By comparing the coefficients of the Fj's in ([30} and d26 



( 27 1 one finally obtains a differential relation between the 
gi's and the controls u\ and u 2 in the form of a matrix 
equation 



1 





sin(2 52 ) 


9i 




U\ 





cos(2<7i) 


-cos(2g 2 ) sin(2gi) 


9i 




u 2 





sin(2 5 i) 


cos(2<7 2 ) cos(2<?i) 


93 








(34) 



the coefficient matrix being the same as in [2], Eq. (3.7). 
The relation g <-> u is invertible if the determinant of the 
coefficient matrix denoted by T> is different from zero 



Then we have 



cos(25 2 ) 





\V\ = cos(2 ff2 ) ^ 0. 



v -i 



cos(2# 2 



(35) 



(36) 



-sin(2gi) sin(2g 2 ) cos(2gi) sin(2.g 2 ) 
cos(2gi) cos(2g 2 ) sin(2pi) cos(2g 2 ) 
— sin(2pi) cos(2gi) 



g = V^u, 



where u and g are defined by 



(37) 



(38) 



IV. MODEL PARAMETRIZATION 

Because the system has two (scalar) controls we can 
choose two of the three base functions gi freely correspond- 
ing to free selection of the two controls. The third base 
function has to be determined from the last equation of 
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u 2 


, 9 = 


92 









93 



(34 1. Parametrization actually means that the input-output 
pairs can be determined from the parameter functions without 
explicitly solving of the system equations according to Fig. 
1. Due to the flatness-based design idea, computation of 
the third base function as well as of the controls must not 
include integrations as given by the equations Q. Only 
differentiations are allowed. Consequently, based on the third 



equation in (34 1, the base functions g 2 and 53 are chosen 
as parameter functions. Then these are also so-called fiat 
outputs, see [8], denoted by z — (zi,z 2 ) — (^,53). The 



Input u(t) 



Output y(t) 




Parameter function z(t) 

Fig. 1. Mappings 1 & 2 give the corresponding input-output pairs (u,y) 
without explicitly solving the system equations. 



parametrization obtained in this way for g\ and the controls 
are given by 



1 

gi = — arctan 



cos(2g 2 ) ^ 

92 



(39) 
(40) 

9 3 2 cos2 (2. 92 ). (41) 
The state variables are calculated by using the state transition 



ui = .91 +.9 3 sm(2sr 2 ) 

«2 = 



matrix equations ( 23 1 and ( 24 1 



x(t) = Q)x(0) = e 9ltl e 92t2 e S3i * 3 x(0). 



gffa-Fb e S3^3 

V. Control objective 



In population transfer problems from the level 1 corre- 
sponding to the situation 



\M0)\ 2 = Mo) 2 + X2 (o) 2 = 

to the level 2, where 

|V>2(T)| =x 3 (T) 2 +x 4 (T) 2 = 0, 



(42) 



(43) 



where T is the transfer time, we can parametrize the partial 
trajectory by using a sufficiently smooth, but otherwise 
arbitrarily chosen, parametrization xi,x 2 with the boundary 
conditions 



Zl(0) 2 +Z 2 (0) 2 : 

Xl (T) 2 + x 2 (T) 2 



0, 



By dividing the state vector into two parts 

x(t) = (w(t),v(t)) 

w(t) = ( Xl {t),X 2 (t)) 

v(t) = (x 3 (t),x 4 (t)) 



(44) 
(45) 



(46) 
(47) 
(48) 



we can represent the task of driving the state from the initial 
one to the final one in a finite time T as follows 



x(0)= " x f =x{T) m 



We must have D = due to the requirement Dv n = 
for arbitrary vq = ( £30 , £40 ) satisfying the requirement 
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cos/3 
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sin /3 


sin a 







cos a 









1. Then we have two alternatives in (59): 




L»w = 0. (62) 



(50) 



We have chosen a specific parametrization for the initial 
and final values of the state, because the sum of the squares 
of the nonzero state components must be equal to 1 at the 
both ends of the planned trajectory. 

VI. PARAMETRIZATION DESIGN 

The state transition equation x(T) = Q(T, 0)x(0) can now 
be written in the form 



ci = s 2 

or D = 

b) s\ = c 2 = 

These conditions are obtained from the two basic alternatives 

7T 



lcoBfli(T)=0, gi(T)= 2 
\sin ff2 (T) = 0, 52 (T)=0, 



ft) 



'mng 1 {T) = 0, 9l (T) 
cos ff2 (T) = 0, 52 (T) 





7T 



In the case of the first alternative a) we have 

\Si = sin gi (T) = 1 
lea =cos ff2 (T) = 1. 



(63) 



(64) 



(65) 



Wt 




' A 


B ' 
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D 







(51) Consequently, 



wt = Bvq 
= Dv . 



D 



(52) 



where A, B, C, and D are 2 x 2-blocks of the 4x4- 
dimensional state transition matrix $(T, 0). 
For the state transition matrix 



wt = Bvq = 



-S1C2 



sin 53 
cos 33 



.S3 
C3 



-C3 
S3 

cos 53 
sin 53 



sin a 

cos a 



$(t,0) = e 9lFl e 92i?2 e 93F3 



(53) 



sin 53 sin a — cos 53 cos a 
cos 53 sin a + sin 53 cos a 



^ ft 

fc=0 



we obtain the series representations in closed form 



cos gi I + sin 5^ Fi 



(55) 



(66) 

(67) 
(68) 
(69) 
(70) 

In the same way the alternative b) can be solved giving 

9z(T) = l -(a + 0). (71) 



series 






cos(-.g 3 - a) 
sin(-5 3 - a) 




cos j3 
sin j3 
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(54) 




■ : 9s(T) = - 


(q 4 





due to the fact that Ff = —J, i = 1,2,3, where I is 4x 
4 identity matrix. Then the product of the three exponent 
functions is of the form 

$ = ( Cl J + SiFiXcaJ + s 2 F 2 )(c 3 / + s 3J F 3 ) (56) 
c,;=cos5i, s,; = sin^i, i= 1,2,3. (57) 

Now the D-part and S-part of the transfer matrix $ are given 
by 

D = 



Due to trigonometric functions in the equations there are also 
other possibilities for the final values of g 2 and 53 deviating 
by the multiples of it or 2ir. These possibilities need further 
considerations and are not studied here. We choose the 
alternative b) for the basis of our control design. So, we 
have to find sufficiently differentiable parameter functions 
g 2 and 53, which together with the dependent basis function 
<?i have to satisfy the boundary conditions 
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= CiS 2 
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S3 
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S3 


• (61) 



5i(0)=0, 5 i(T) = 
52(0) =0, 52 (T) = ^ 
33(0) =0, 53 (T)=" 



(72) 



(a + 0). 



The final value of g\ depends on the derivatives of g 2 and 
53. This means that we have to adjust these derivatives via 
the equation (39 1 to agree with the requirement gi{T) = 0. 

Carefully planned and realized simulations are needed to 
confirm the feasibility of our parametrization approach. 



VII. SIMULATION STUDY 

First preliminary simulation results demonstrate that the 
methodology developed actually drives the state of the 
system from the given initial state (level 1) to the given 
final state (level 2). A minimal parametrization for the 
parameter functions 52 and 173 were chosen without any 
specific optimization procedure. The only requirements are 
that the given boundary conditions f72| ) are satisfied, and that 
the equation which gives the base function g\ also gives the 
correct initial and final values for g\. The following values 
were used in the simulations 



a =-2tt/3 : 
(3 = vr/3 , 
T =10. 



Then the final value for 53 becomes 

•57T 



03 cn 



2.62. 



(73) 



(74) 



Because g% has to change from to n/2, we chose the linear 
function 



92{t) = 
The boundary values 



7T t 

2 T 




0.157*. 



= 
= 

are obtained when we choose 

= 
= ' 

Then the third order polynomial suffices 

t / t \ 2 




t \ 3 



93(t) = ao + Ol j, + ^{jf) +fl3 (y) 



(75) 



(76) 



(77) 



(78) 



The coefficients are obtained from the boundary conditions, 
giving finally 



*»(*)= 7 |3(j;) 
The binding condition 



/ \ 2 / t \3 

-2 - 
T 



7 = 2 - (« + P) ■ ( 79 > 



92 sin(2ffi) + g 3 cos(25 2 ) cos(2#i) = (80) 
gives the base function g\ for the given parameter functions 



(75 1 and (79 1. The functions are depicted in Fig. 2 and 3. 



The controls were calculated by using the formulas (40 1 and 



(41 1. They are depicted in Fig. 4. 



The behaviour of the state variables are given in Fig. 5 
indicating that the desired final state, where 

x 3 (T) =Xi(T) =0 

has been obtained. The simulations were carried out and the 
figures produced by using Mathematica 7 package [25]. 






Fig 



2468 10 02468 10 
2. The independent parameter functions: 92 - left. 93 - right. 
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Fig. 3. The dependent base function g\. 




Fig. 4. The control variables u\(t) and U2(t) for t E [0, 10]. 




Fig. 5. The state variables obtained by simulating the system equations 
by using the given control functions. Left: x± - lower, X2 - upper. Right: 
x-j - lower, X4 - upper. 




Fig. 6. Population change: x\ + x\ increasing, x| 



VIII. CONCLUSIONS 



x\ decreasing. 



The parametrization idea for constructing open-loop con- 
trols for uncontrollable bilinear systems is applied here. 
We have also earlier studied parametrization of systems 
described by partial differential equations and pseudo- 
differential operator models, see [12]— [14]. Flatness-based 
ideas, originally developed by Michel Fliess and his co- 
workers [7]-[9] have been developed for open-loop control 
design. In some quantum control problems, where laser 
pulses are used for the control, the dynamics is so fast 
that, at least at the present level of the speed of possible 
computations, feedback control seems to be impossible to 
implement even if so-called homodyne detection principles 



can be applied to obtain closed-loop controls. 

Here we studied a two-level population transfer problem. 
Without more advanced differential geometric considera- 
tions, which might be helpful in understanding quantum phe- 
nomena in general, we use the formulation found generally 
in the literature, to obtain our basic driftless system model 
of the form x = g(x)u, where g is linear in the state x. 

Simulation study was required to confirm the quantum 
control approach chosen. Then depending on the choise of 
the alternatives a) or b) different state trajectories can be ob- 
tained resulting, however, the same final state of the system 
when the flatness-based control is applied. Our preliminary 
simulations were based of the alternative b). 

The basic technique applied here is useful also in multi- 
qubit systems and in controlling entanglement of, say, two or 
more qubits. Then tensor product formalism in the Euclidian 
framework is a feasible alternative in the system model 
design. 
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